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Abstract 

The meaning of the inflationary slow-roll approximation is formalised. Comparisons 
are made between an approach based on the Hamilton- Jacobi equations, governing 
the evolution of the Hubble parameter, and the usual scenario based on the evolution 
of the potential energy density. The vital role of the inflationary attractor solution 
is emphasised, and some of its properties described. We propose a new measure of 
inflation, based upon contraction of the comoving Hubble length as opposed to the 
usual e-foldings of physical expansion, and derive relevant formulae. We introduce 
an infinite hierarchy of slow-roll parameters, and show that only a finite number of 
them are required to produce results to a given order. The extension of the slow-roll 
approximation into an analytic slow-roll expansion, converging on the exact solution, 
is provided. Its role in calculations of infiationary dynamics is discussed. We explore 
rational-approximants as a method of extending the range of convergence of the slow- roll 
expansion up to, and beyond, the end of infiation. 
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1 Introduction 



Inflationary universe models are based upon the possibility of slow evolution of some scalar 
field 6 in a potential y((/>) Although some exact solutions of this problem exist, most 

detailed studies of inflation have been made using numerical integration, or by employing 
an approximation scheme. The 'slow-roll approximation' [P, |5|, which neglects the most 
slowly changing terms in the equations of motion, is the most widely used. Although this 
approximation works well in many cases, we know that it must eventually fail if inflation 
is to end. Moreover, even weak violations of it can result in significant deviations from 
the standard predications for observables such as the spectrum of density perturbations or 
the density of gravitational waves in the universe |^, As observational data sharpen, 
it is important to derive a suite of predictions for the observables that are as accurate as 
possible, and which cover all possible inflationary models. 

In the literature, one finds two different versions of the slow-roll approximation. The 
first 1^, ^ places restrictions on the form of the potential, and requires the evolution of 
the scalar field to have reached its asymptotic form. This approach is most appropriate 
when studying inflation in a specific potential. We shall call it the Potential Slow-Roll 
Approximation, or PSRA. The other form of the approximation places conditions on the 
evolution of the Hubble parameter during inflation 0]. We call this the Hubble Slow-Roll 
Approximation, or HSRA. It has distinct advantages over the PSRA, possessing a clearer 
geometrical interpretation and more convenient analytic properties. These make it best 
suited for general studies, where the potential is not specified. 

In this paper, we clarify the meaning of the different slow-roll approximations that exist 
in the literature, which often describe a variety of slightly different approximation schemes 
applied to different variables at different orders. By formalising the slow-roll approximation 
in detail, we will show how to use it as the basis of a slow-roll expansion — a sequence of 
analytic approximations which converge to the exact solution of the equations of motion for 
an inflationary universe. Such a technique relies strongly on the notion of the inflationary 
attractor, whose properties we describe. The use of Fade and Canterbury approximants 
^ allows us to further improve the range and rate of convergence of this slow-roll expansion. 



2 Equations of Motion and their Solution 

We shall deal with the equations of motion in two different forms, both appropriate for a 
homogeneous scalar fleld 4>, evolving in a potential V{4')- We assume that enough inflation 
has occurred to render the densities of all other types of matter negligible, and to establish 
homogeneity in a patch at least as big as the horizon. The most familiar form of the 
equations, in a zero-curvature Priedmann universe, is 

4> + 3H,j} = -V , (2.2) 

where H = a/a is the Hubble parameter, a is the scale factor (synchronous rather than 
conformal), mpi the Planck mass, overdots indicate derivatives with respect to cosmic time 
t, and primes indicate derivatives with respect to the scalar field (j). 
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One can derive a very useful alternative form of these equations by using the scalar field 
as a time variable |l^. This requires that does not change sign during inflation. 

Without loss of generality, we can choose > throughout (this will determine some signs 
in later equations). Differentiating Eq. (|2.1|) with respect to t and using Eq. (|2.2D gives 



2H = --^^^. (2.3) 



m 



PI 

We may divide each side by (j) to eliminate the time dependence in the Friedmann equation, 
obtaining 

[H'm'-^H\<t>) = -^^V{^), (2.4) 



m 



PI 



H'icP), (2.5) 
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and Eq. (p.3|) implies H < 0. This new set of equations — the Hamilton-Jacobi equations 



— is normally more convenient than the Eqs. ( ^.1| ) and (2^). They were used by Salopek 
and Bond Q to establish several important results to which we refer later. 

These equations allow one to generate an endless collection of exact inflationary solutions 



via the following procedure ||T^, 11]: choose a form of H{(f)), and use Eq. (p.4[) to find 



the potential for which the exact solution applies; now use Eq. ( p. 51) to find which 
allows the (/)-dependences to be converted into time-dependences, to get H{t); if desired, 
a further integration gives a{t) (though this last step is seldom required). For example, 
this procedure gives a very easy derivation of the exact solution describing 'intermediate' 
inflation, which corresponds to the choice H{(j)) oc with a a positive constant and 
> 0. This exact solution was derived using the Hamilton-Jacobi equations by Muslimov 
]I0| (and independently by Barrow ||T^ using a different technique). Other fully integrated 



exact solutions have been found by Barrow |14, 15| 



It is normally impossible to make analytic progress by first choosing a potential V{4')^ 
because Eq. (^^) is unpleasantly nonlinear. The simplest exception is the exponential 
potential, known to drive power-law inflation, for which the Hamilton-Jacobi formalism was 
used by Salopek and Bond 0| to find (in parametric form) the general isotropic solution. 

2.1 The potential-slow-roll approximation (PSRA) 

When provided with a potential V{(t)) from which to construct an inflationary model, the 
slow-roll approximation is normally advertised as requiring the smallness of the two param- 
eters (both functions of ^), defined by 



IGtt V V{(j)) ) 



PI 



(2.6) 
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Henceforth, we refer to them as potential- slow-roll (PSR) parameterJ^ Their smallness is 
used to justify the neglect of the kinetic term in the Friedmann equation, Eq. and 
the acceleration term in the scalar wave equation, Eq. ( |2.2| ). Unfortunately, the smallness 
of the PSR parameters is a necessary consistency condition, but not a sufficient one to 
guarantee that those terms can be neglected. The PSR parameters only restrict the form 
of the potential, not the properties of dynamic solutions. The solutions are more general 
because they possess a freely specifiable parameter, the value of 4>, which governs the size 
of the kinetic term. The kinetic term could, therefore, be as large as one wants, regardless 
of the smallness, or otherwise, of these PSR parameters. 

In general, this PSR formalism requires a further 'assumption'; that the scalar field 
evolves to approach an asymptotic attractor solution, determined by 

The word 'assumption' is placed in quotes here because, in general, one is able to test 
whether Eq. (2.8) is approached for a wide range of initial conditions. The inflationary 



attractor is of vital importance in the application of the slow-roll approximation, and we 
discuss its properties in Subsection |2.3|. 



2.2 The Hubble-slow-roll approximation (HSRA) 

If H{(j)) is taken as the primary quantity, then there is a better choice of slow-roll parameters. 
We define the HSR parameters, e„ and r/^, by 



These possess an extremely useful set of properties which make them superior choices to €v 
and rjY as descriptors of inflation: 



^TT^^ =—n-rU (2-11) 



We have exactly 



y + 02/2 V dlna 

(f> (_ dln(j) _ dlnH'\ 

ejj <C 1 is the condition for neglecting the first term of eq. ( |2.4| ) [the kinetic term in 



Eq. dp)] 



^To preview what is to come, these parameters are sufficient to obtain results to first order in slow-roll. 
However, the general slow-roll expansion requires an infinite hierarchy of parameters which will be defined 
in Section ^ 
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• |r/j^| <C 1 is the condition for neglecting the derivative of the first term of Eq. ( p. 41) 
[the acceleration term in Eq. (|2.2| )]. As a consequence, all the necessary dynamical 
information is encoded in the HSR parameters. They do not need to be supplemented 
by any assumptions about the inflationary attractor, Eq. ( |2.8D . 

• The condition for inflation to occur is precisely 

a>0^^ejj<l. (2.13) 



There is an algebraic expression relating to en and ij^ (using Eq. (|2.4|)): 



The true endpoint of inflation, gauged by the HSR parameters, occurs exactly at = 1. 
When using the PSR parameters, this condition is approximate; inflation ending at 6^ = 1 
is only a first-order result. 

For rjy, the relation to the HSR parameters is differential rather than algebraic, 



although a more compact representation in terms of higher-order parameters will be pre- 
sented in Section ^. This will show that the first term in Eq. ( 2.15| ) is of higher-order in 



slow-roll, so that to lowest-order, one has rjy = tj^ + ejj. Note that 7]jj and rjy are not the 
same to first-order in slow-roll, as one expects from oc V. We could have defined r]H 
to coincide with rjy in slow-roll, by defining fj^ = rjff — e^, but we prefer to regard the 
definitions in Eqs. ( p.9[ ) and ( 2.10| ) as fundamental. 



The definitions can be used to derive two useful relations between parameters of the 
same type 

Note that although, as functions, the parameters of a given type, either HSR or PSR, are 
related, their values at a given (f) are independent of one another. One immediately sees the 
different 'normalisation' of the r] from Eqs. ( p. 16 ) and ( 2.17| ). 



It is important to stress that although we have derived self-consistent exact expressions 
relating the PSR to the HSR parameters, we cannot invert these expressions without first 
assuming that the evolution has reached the attractor, Eq. (|2.8D . As already mentioned, 
the attractor constraint is part of the structure of the HSRA, but is absent from the PSRA. 
So, while the HSRA implies the PSRA, the converse does not hold without assuming the 
attractor constraint. 
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2.3 The inflationary attractor 

Already in this Section we have seen how important the notion of the inflationary attractor 
is. The behaviour of this attractor was established by Salopek and Bond Q, and we now 
discuss its properties. 

Suppose HQ{(j)) is any solution to the full equation of motion, Eq. (|2.4D — either in- 
flationary or non-inflationary. Consider, first, a linear perturbation 5H{(j)). We shall also 
assume, and discuss further below, that the perturbation does not reverse the sign of 0. It 
therefore obeys the linearised equation 

127r 

HL5H' ~ -^Hq5H , (2.18) 

^Pi 

which has the general solution 

5H{4>) = 5H{(t>,) exp ^64] . (2.19) 

ympi J(t)i J^o J 

Since Hq and dcp have, by construction, opposing signs, the integrand within the exponential 
term is negative definite, and hence all linear perturbations die away. 

If Hq is inflationary, the behaviour is particularly dramatic because the condition for 
inflation bounds the integrand away from zero. Consequently one obtains 

<5F(0) < 6H{<Pi) exp ( 10 - 1 . (2.20) 

That is, if there is an inflationary solution all linear perturbations approach it at least 
exponentially fast as the scalar field rolls. 

Another way of writing the solution for the perturbation, regardless of whether Hq is 
inflationary or not, is in terms of the amount of expansion which occurs, by using the 
number of e-foldings N as defined in the following section (Eq. ( |3.lD ). We get the precise 
result g 

5H{(j)) = 5H{(t)o) exp {-^Ni - N\) . (2.21) 

For non-linear perturbations, the problem is more complex; though all the solutions 
are easily seen to approach each other, we have not shown that they do so exponentially 
quickly. The most awkward case is where a perturbation actually reverses the sign of as 
the Hamilton-Jacobi equations are singular when that happens. Nevertheless, as long as the 
perturbation is insufficient to knock the scalar field over a maximum in the potential, the 
perturbed solution will inevitably reverse and subsequently pass through the initial value 
again; then it can be treated as a perturbation with the same sign of (j) as the original 
solution. 

The picture that emerges is therefore as follows. Provided the potential is able to sup- 
port inflation, the inflationary solutions all rapidly approach one another, with exponential 
rapidity once in the linear regime. Even when inflation ends, the universe continues to 
expand and therefore the solutions continue to approach one another. Consequently, even 
the exit from inflation is independent of initial conditions. Note that there is no concept 
of a single 'attractor solution': all solutions are attractors for one another and converge 
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asymptotically. As we shall see, this is a vital requirement if a slow-roll expansion is to 
make any sense. 

The situation where the inflationary attractor does not apply is therefore soon after 
inflation begins. Normally this is in the distant past and of no concern. An exception, 
recently noted is hybrid inflation. There, the slow-roll parameters rise above unity, 
halting inflation, and then fall back below unity, reaching very small values. Nevertheless, 
it is easy to show that inflation fails to restart despite the smallness of the PSR parameters, 
because there is insufficient time for the solution to approach the inflationary attractor. 
Another similar situation will be discussed in Section ^. 



3 A Better Measure of Inflation 

An important quantity for making inflationary predictions is the amount of inflation that 
has taken place. Inflation is commonly characterised by the number of e-foldings of physical 
expansion that occur, as given by the natural logarithm of the ratio of the flnal scale factor 
to the initial one. This can be expressed exactly as. 



N^ln^-l = J^['''d^. (3.1) 

If one is working in the PSRA then , provided the attractor solution Eq. ( |2.8D is attained, 
this may be approximated by 



/ 47r f^f 1 

These formulae are needed to make the connection between horizon-crossing times in cal- 
culations of the production of scalar and tensor perturbations. A comoving scale k crosses 
outside the Hubble radius at a time which is N(k) e-foldings from the end of inflation, 
where 

A^(^)=62-ln4--lni^^^+ln%-iln%. (3.3) 

The subscript '0' indicates present values; subscript 'k' specifles the value when the wave 
number k crosses the Hubble radius during inflation (fc = aH); subscript 'end' specifles the 
value at the end of inflation; and pi-eh is the energy density of the universe after reheat- 
ing to the standard hot big bang evolution. This calculation assumes that instantaneous 
transitions occur between regimes, and that during reheating the universe behaves as if 
matter-dominated. Ordinarily, it is taken as a perfectly good approximation that the co- 
moving scale presently equal to the Hubble radius crossed outside the Hubble radius 60 
e-foldings from the end of inflation, with all other scales relevant to large-scale structure 
studies following within the next few e-foldings. 

Something like 70 e-foldings is normally advertised as the minimum for inflation to 
solve the various cosmological conundrums, such as the flatness and horizon problems. 
However, it is well known that this is an approximation (for example in the standard big 
bang model, the universe expands by more than a factor e^^ between the Planck time and 
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the present, without solving the flatness or horizon problems), based on the assumption of a 
constant Hubble parameter. A better measure of inflation is the reduction in the size of the 
comoving Hubble length, 1/aH. First of all, the condition for inflation (a > 0) is equivalent 
to d[{aH)~^]/dt < 0. Secondly, it is the reduction of 1/aH, not that of 1/a, which solves 
the flatness and horizon problems. And finally, for the generation of perturbations, it is the 
relation of the comoving wavenumber k to aH that is important. We therefore define 

iV^lni^. (3.4) 



It can be shown that 



/ 47r r't'f 1 

Since H always decreases, iV < by definition, with the difference indicating the extra 
amount of expansion, required by the decrease of H during inflation. In the extreme slow- 
roll limit (HSR or PSR) A^ and N coincide. This tells us that the true condition for sufficient 
inflation should be that A" (not A) exceeds 70. 

Again, using the potential, we can only write down an approximate relation. At lowest- 
order, N and A^ coincide, while to next order one obtains 



/ 47r /■<?/ 1 / 1 1 \ 

N{4>^,<|)f)^-J^ / ^== l--6v.(0)--w(</') (3.6) 

Eq. ( ^.5| ) gives A during an arbitrary inflationary epoch (not just quasi-de Sitter). It 
also holds if inflation is interrupted for a period, while the dynamics are still dominated by 



the scalar field. This requires > 1, causing the integrand in Eq. (3^) to change sign over 



a range of (p. This is unlikely, but it has recently been noted |17| that it arises in a variant of 
hybrid inflation with a quartic potential, where the potential that drives inflation is of the 
form V (X {1 + A0^). For some values of A, the potential temporarily steepens sufficiently to 
suspend inflation, while 1/aH increases. This must be compensated by extra inflationary 
expansion later on. 

To make use of the new expressions, Eqs. ( p.5[ ) and ( |3.6| ), we need an A(A;) relation to 
replace Eq. (|3.3D . This is simply 

m = 62 - In - In 1^ " ^ 1- ^ • (3-7) 



Although this strongly resembles Eq. ( |3.3D , it is in fact slightly simpler; the difference 
between Vfc and 14nd is now part of the definition of A. 



4 The Hierarchy of Slow-Roll Parameters 

We now reconsider the intrinsic structure of the slow-roll approximation. The HSR pa- 
rameters (e^ and rjfj of Section p. 2D , which measure the first and second derivatives of the 
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Hubble parameter, are all we require to obtain results to first-order in the slow-roll expan- 
sion. However, to go beyond this, we require more derivatives, necessitating further slow-roll 
parameters. In general, there will be an infinite number of these, incorporating derivatives 
to all orders; we shall prove that each additional order in the expansion requires the intro- 
duction of one new parameter. The formal order of the expansion parameter depends on 
the number of derivatives it contains, two derivatives for each order. 

In [^], where second-order results were derived, the extra parameter was defined as 
^CKLL = {^pi/^''^)H"' / H' , which is of the same order as the others. This definition is rather 
unfortunate; we are meant to be expanding about a flat potential, and this parameter is 
not guaranteed to tend to zero in that limit because of the derivative in the denominator. 
This definition has led to some confusion (for instance, in ||l^ where ^ was not treated as an 
expansion parameter at all). The moral of this is that the hierarchy of slow-roll parameters 
should be carefully defined, so as to tend to zero as the potential approaches flatness in 
arbitrary ways. Even with this restriction, there are different ways one could define the 
hierarchy. If a superscript indicates the n-th derivative with respect to (p, the simplest 
definition would appear to be 



m 



PI 



H 



2/n 



(4.1) 



which gives ^/Sj, 



€h and 



rjjj. However, there is a superior alternative. 



4.1 The HSR hierarchy 

We shall work with the set of definitions 




which can be expressed as 



din a 



(4.2) 



m 



(4.3) 



These have elegant properties, which we shall use to recast the HSR approximation as 
an HSR expansion. Furthermore, as we demonstrate below, only a finite number of these 
parameters are needed to obtain results to any given order. 

It is difficult to incorporate into this scheme naturally — it must be defined separately. 
We use the form given in Eq. (p.9|), but shall refer to en as ^j3H in later Sections. The first 
parameter this definition yields is ^{3^ = Vh, in accord with Eq. ( ^.10 ). The next four HSR 
expansion parameters are 



m 



PI 



fH'H'' 



4tt V 

mil fH'^H' 



(4.4) 
(4.5) 
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ru^'Pu = ^(^^^1 , (4.6) 



mil 
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mil 
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The strength of these definitions is that each parameter combines H and its derivatives, 
raised as a whole to some power 1/n, where n E Z~^. This has the interesting consequence 
that if we wish to convert the Taylor-series of some function of H{(p), and its derivatives, 
into an expansion in slow-roll parameters, we are guaranteed that the powers of any specific 
parameter, "/Jj^, will be integer multiples of n. Thus, the lowest-order term we can expect to 
find involving this parameter will be {"'Ph)^- Consequently, if one is interested in expanding 
this function up to order m, then at most the first m + 1 parameters can appear. Thus, 
despite the potentially infinite number of slow-roll parameters, we require only a finite 
selection of them at any order. 

4.2 The PSR hierarchy 

In the spirit of the previous subsection, we introduce a hierarchy of PSR parameters. As 
stressed earlier, these are not as useful in calculations of infiationary dynamics as the HSR 
parameters. They only classify the fiatness of the potential, and so encode no informa- 
tion about initial conditions. This necessitates adding the additional attractor constraint, 
Eq. (pj). 



As before, we retain the standard definitions of and rjy from Eqs. (2^) and ( p. 71) , 
and encapsulate the definition of r]y, and all higher-order parameters, in the new set of 
quantities 



i=l 



(4.8) 



which reduce to 



m 



2 



This allows construction of a set of PSR expansion parameters 

^v= Pv = ^(^— ' (4-11) 
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In Eqs. (2.14)-(2.17), we presented some additional properties of e and rj. As one 



may suspect, these can be written in terms of the higher-order parameters. An extensive 
collection are given in the appendix; here, we note the exact relations 

ev = end^y , (4.14) 



3 - e 



H 



W = + . (4.15) 

3 — fin 

The latter demonstrates that rjv = €h + tjh in the lowest-order HSRA, as stated earlier. 



5 From Slow-Roll Approximation to Slow-Roll Expansion 

A procedure for generating analytic approximations to inflationary solutions, from a slow- 
roll approximation, has a broad spectrum of applications. More useful still is a slow-roll 
expansion, from which solutions could be generated analytically to any required order in 
the slow-roll approximation. An indication of how to go about this was given by Salopek 
and Bond |Q. We now show how this may be achieved within the framework of the slow-roll 
parameters. 

Let us first emphasise the importance of the attractor behaviour for this procedure. The 
general isotropic solution for a given potential possesses one free parameter, corresponding 
to the freedom to specify H (or equivalently (p) at some initial time. However, the traditional 
slow-roll solution, and its order-by-order corrections that we shall describe, generate only 
a single solution. Unless an attractor exists and has been attained, there is no need for 
this single solution to represent in any way the true solution for that potential. However, 
if the attractor has indeed been reached, then any particular solution provides an excellent 
approximation to those arising from a wide range of initial conditions. This is particularly 
important when inflation approaches its end, and the slow-roll parameters becoming large, 
because one might naively assume that the one-parameter freedom could be important 
there. If the attractor solution exists, then solutions for a wide range of initial conditions 
will converge, and subsequently all exit inflation in the same way. Hence, an expansion 
approximating one particular solution serves as an excellent approximation to them all, 
provided the initial condition for that solution is not pathological (which is prevented by 
the assumption that energies are less than the Planck energy). 

We note that there is a formal problem in attempting to prove 'no-hair' theorems; if 



inflation is to end there is no formal asymptotic regime [14|. However, our requirement is 
just that enough inflation occurs to ensure that the range of values of (j) needed to encompass 
the entire family of solutions is sufficiently small to validate the Eq. ( |2.8D . In typical models 
so much inflation occurs that this situation is easily achieved. 

5.1 The traditional Taylor-series approach 

We start with a potential V{(j)) for which a solution is desired. Typically, we cannot solve 
exactly for H{(p); instead, we aim to find an approximate solution, in terms of V{(f)) and 
a multivariate Taylor expansion in the PSR parameters, all of which can be computed 
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analytically given an analytic V{(j)). First recast Eq. (2.4) as 



Svr 
3mpj 



V{4>) 



1 - 3.„( 



Then seek an approximate solution for H'^{cl)) of the form 
Svr 



3m% 



y((/)) (l + aev + brjv + eel + devr]v + er?^ + /'^v + 



(5.1) 



(5.2) 



where a, 6,... are constants to be determined. [In fact, Eq. ( 2.14 ) already guarantees that 
6 = e = / = 0, by forcing every term in to contain at least one power of e^.] Note 
also that, for reasons discussed in Section 4.1, is the lowest-order term involving to 
appear in Eq. (|5.2|) , and that no further PSR parameters appear at second-order. Thus, 



we require an expansion of (1 — 6^/3)" in PSR parameters of ascending order. This is 
most readily achieved by assuming general expressions for the HSR parameters, in terms of 
PSR parameters, with a similar form to Eq. ( |5.2| ). Starting with general first-order forms, 
these may be substituted into Eq. ( p. 14 ) , which is then solved for the unknown constants by 
comparing coefficients|^. After first-order results have been obtained, this procedure may be 
repeated iteratively order-by-order. It is a straightforward (albeit tedious) matter to invert 
the expansion of (1 — e^^/S) using the binomial theorem, and hence obtain the corresponding 
series for i?^(0). 

We have done this to fourth-order, although the results may be truncated at lower-order 
to obtain more manageable expressions, as desired. We find 



3 

62 2 



. . 14 3 



32 o 5 



10 



16 



.2 a2 



9 '^^^ + 3' 



-vVv + 

, 2 
^+27^ 



-Ml 



44 



+ 



530 



9 " 27 

The second-order truncation concurs with the result derived in |19| 
Hence, we find the approximate solution 



+ O5 



(5.3) 



Svr 

2,171% 



327 



1 



25 
— ( 
27 



+ 



27 



SI 



+ 



460 

"sT 



172 



■2 J2 



14 



SI ^ SI 



26 
~ 27 
44 
81 



g2 .2 



(5.4) 



Thus, we have generated an analytic solution for inflation in the potential ^((/>), that is 
accurate up to fourth-order in the slow-roll parameters, rather than the usual lowest-order. 

We illustrate this with the specific example of the potential V{<j)) = m^(/>^/2. To second- 
order in the HSR parameters, one finds 



m 



PI 



-H 



4vr(/)2 



1 



m 



PI 



+ 



m 



PI 



(5.5) 



^Although the slow-roll functions are all inter-dependent, as all are based on derivatives of V{(fi), their 
values at a given (j> are independent. Consequently, one should imagine that this procedure is being carried 
out separately at each (j). The results are then used to construct a function of cj>. 
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and hence 

^<*>=^ 

Results accurate to one order less, as given in Ref. can be obtained by removing the 
last term from both expressions. The behaviour of H{(j)) is indicated in Figure 1, where for 
comparison, the exact numerical solution is also shown. 

For this particular potential, it is necessary to include both first and second-order cor- 
rections in order to maintain, as a sensible definition of the end of inflation, the condition 
ejj{4>) = 1 (here, we refer to the precise €^{4') derived from the approximate solution H{(l)), 
and not to enifp) truncating to the same order). To lowest-order, this is guaranteed because 
the potential has a minimum where V{(f)) = 0; but, because corrections become large near 
the end of inflation, they may spoil this. In fact, for a (/)^-potential, if one includes only the 
first-order corrections, they conspire so that 6^(0) fails to reach unity for any </), despite the 
solutions being closer to the exact numerical solution than the lowest-order results for the 
bulk of the evolution. Including the second-order corrections removes this problem for the 
0^-potential. The exact behaviour of e^((/)) is shown in Figure 2 in each case. For potentials 
with (j)" behaviour with a > 4, the end of inflation is well defined, even at first-order. 

5.2 The rational-approximant approach 

In the chaotic inflationary example presented above, it was seen how, as the field rolls 
toward the minimum of the potential, the Taylor-series expansion becomes progressively 
less accurate due to the dependence of the slow-roll parameters on inverse powers of (p. It is 
not true, in general, that eH{4>) will be the first parameter to become large in this manner. 
Hence, it is not guaranteed that the approximation will describe the evolution through, or 
even up to, the end of inflation. At (j) = 0, where V{4>) = 0, the Taylor expansion diverges, 
although inflation necessarily finishes before this. 

In cases such as this, rational-approximant techniques |8|, ^ can be effective. For the 
single variable case, instead of using a single Taylor-polynomial, we approximate using the 
Fade approximant — a quotient of two polynomials — in the hope of achieving a better 
range and rate of convergence. A particular application of this technique to infiation was 



made in Ref. |20|. Unfortunately, Fade approximants can only be directly applied to single 



variable expansions, and we require the extension of this theory to multi-variable problems. 
5.2.1 The Canterbury approximant 

The Canterbury approximant Q supplements a Fade quotient approximant in many vari- 
ables with a minimal Taylor series — where minimal means containing as few terms as 
possible. So, for a function of r variables f{xi,X2, ■ ■ ■ ,Xr) expanded to n'^^-order, we pos- 
tulate an approximant of the form 



A / \ n n n r 

^ ^' ^' ' '/ pj=0p2=0 Pr=0 Ll=l 



(5.7) 



where A{xi,X2, ■ ■ ■ , x,.) and B{xi,X2, ■ ■ ■ , a^r) are multi-variable polynomials of order L and 
M respectively, the ep^pj...^^ are constants and pi are the powers of Xj, ranging between and 
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n. The orders of the approximants, L and M, are constrained by the relation L + M = n, 
and the polynomial B{xi,X2, ■ ■ ■ ,Xr) contains no zeroth-order term. The undetermined 
constants in the approximant are fixed by expanding [L/M]f to n^^-oidev, and matching it 
with the n*^-order expansion of /(x). Such an approximant often yields faster convergence 
because the quotient contains better estimates of the higher-order terms than a truncated 
Taylor series. Ref. provides a useful introduction to approximant techniques. 

This approximant is not unique. Any attempt to neglect the Taylor terms in Eq. ( ^.7[ ) 
and solve for the remaining constants fails in general, as unfortunately there will be an 
insufficient number of free constants to do this. The choice of which values we assign to 
the constants in the Pade-term (and hence the form of the corrective Taylor series), is in 
fact completely arbitrary, demonstrating the non- uniqueness of Eq. ( |5.7D . We find however, 
that keeping the Taylor-correction terms purely n'^^-order simplifies the analysis somewhat, 
and also ensures that they are small in slow-roll. 

We are now in a position to recast the PSRA in the form of Eq. (|5.7| ). In terms 
of PSR parameters, during slow-rolling inflation, an arbitrary function f{H{(j))) may be 
approximated by 



A ( La \ n n n 



POPl---Pn 



po=Opi=0 p,i=0 



i=0 



(5.8) 



where the "Pv are the general PSR parameters defined in Eq. (|4.8|); recall, also, that we 
have ^Pv = (-V- 

This complicated formalism is clarified by showing it at work to a given order. We 
present a [2/2] Canterbury approximant to (1 — e^/S), and indicate how this may be used 
to construct the approximant to -ff^ (</>). For (1 — 6^/3), we obtain 

1 _L 21 _7 _ 53f: „ _L 89f:2 , „2 _ 2 e2 r, 
[^/^J(1-.h/3) - 1 ,67 _7 _ 7 ,35 2 ,„2_2t2 81 



162 324 ^ 162 27 ^ 



Diagonal Canterbury approximants (ie, the [L/L] cases), share many of the useful prop- 
erties of standard Pade approximants (duality, homographic invariance, unitarity, etc), and 
substantially simplify the application of the Canterbury technique . Use of diagonal ap- 
proximants, where possible, is thus recommended. In particular, the duality property may 
be exploited to save considerable effort when calculating the corresponding expression for 
H'^{4>) from Eq. (^.S|). The duality property is as follows; if / (x) = [g {x)]~^ and g (0) /O, 
then 

We may thus obtain a [2/2]ff2(^^-j from Eq. (|5.9|), via 



imHH^) = ^^('^) {[2/2](i-.W3)}"' • (5-11) 
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5.2.2 Simplified Canterbury approximants 

The Canterbury approximant provides a powerful technique for calculations of inflationary 
dynamics. However, in its full higher-order glory, it can be quite cumbersome and unwieldy. 
It is therefore useful to find circumstances in which the corrective Taylor series is not 
required. The simplest way to bring this about is to take the [0/n] approximant which, 
as we will show, never needs correcting. However, it may also be true that at low orders 
there is sufficient freedom in the diagonal approximants, due to the vanishing of some of 
the terms in the original Taylor series. 

In fact, the [1/1] approximant has this property. Assuming a general form for the 
[l/l]i_^^/3 approximant, and matching to the Taylor expansion for (1 — e/^/S) from Eq. ( ^.3[ ) 
to second-order, yields 

n/Ti _ a + bev + criy ey ( 4 2 . \ 

Comparing coefficients, we arrive at the result 

ivii,.-../3, = i'+ c ^: ■ <"^' 

^ ' 3'/v 

which, by construction, agrees with the Taylor series to second-order. The corresponding 
H{(p) is easily obtained. For the (jy^ potential, examined earlier, this does indeed improve 
on the second-order Taylor series given in Eq. (|5.6| ), as shown in Figures 1 and 2. 

As a final observation, note that if the slow-roll parameters all have the same functional 
form, permitting us to write 

"/3, = ^C,r(<A) (5.14) 

i 

where Cj are constants and /(</)) is an arbitrary function of (f) (which is small in the PSRA), 
then we can always circumvent the need for a corrective Taylor-series by expanding in 
powers of /{(p), thus reducing the problem to the unmodified Fade case. 

6 Conclusions 

By defining a suitable hierarchy of parameters, we have extended the slow-roll approxima- 
tion to a slow-roll expansion, allowing progressively more accurate analytic approximations 
to be constructed via an order-by-order decomposition in terms of slow-roll parameters. 
The use of rational approximants pushes the range of validity of the slow-roll expansion up 
to, and in many instances beyond, the end of inflation. With the accurate observational 
information becoming available, this allows an assessment of the accuracy of calculations 
within the slow-roll approximation, and is especially important with the present consider- 
able emphasis focussed on inflationary models which make predictions far from the standard 
(zeroth-order) case. 

We have used these parameters to define an improved measure of the amount of inflation. 
However, present uncertainties regarding the physics of reheating make it useful only in 
rather extreme circumstances such as a temporary suspension of inflation, during which the 



universe remains scalar field dominated, as in the hybrid inflation model of Ref. [17|. 
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Let us caution the reader regarding the necessity of the attractor condition for the slow- 
roh expansion to make sense. By incorporating order-by-order corrections, we can only 
generate one solution, H{(j)), out of the one-parameter family of actual solutions allowed by 
the freedom of H, or equivalently cj)), permitted by the initial conditions. If the attractor 
hypothesis is not satisfied, then the solution generated — while conceivably an accurate 
particular solution of the equations of motion — need have no relation to the actual dy- 
namical solutions which might be attained. A case in point is the exact 'intermediate' 



inflation solution pO| , 21 1. For small (p, this solution corresponds to the rather unnatural 
(and noninflationary) behaviour of the field moving up the potential and over a maximum, 
beyond which inflation starts. If one attempts to use our procedure to describe this entrance 
to inflation, the solutions generated bear no particular resemblance to the exact solution 
until well into the inflationary regime^. This serves as a cautionary note, that known exact 
solutions are typically only late-time attractors, and unless a significant period of inflation 
occurs before the time of interest, so that the attractor solution is reached, they are of little 
relevance. 

Importantly, with regard to the exit from inflation, we are on much safer ground. It 
is assumed that enough time has passed for the attractor to be reached, and hence all 
solutions exit from inflation in the same way. Therefore, when our expansion procedure 
supplies a particular solution, it provides an excellent description of the way in which the 
entire one parameter family of initial conditions will exit inflation. Without this vital point, 
the generation of solutions via the slow-roll expansion would be fruitless. 

We have concentrated on the dynamics of inflation, rather than on the perturbation 
spectra produced from them. However, the slow-roll expansion can also be brought into 
play there; as an example, we quote the results for the spectral indices n for the density 
perturbations and ut for the gravitational waves (see |^] for precise definitions) . These have 
long been taken as approximately 1 and respectively; results to first-order were given by 
Liddle and Lyth |^] and to second-order by Stewart and Lyth Q. With our definitions, 
these read in the HSRA and PSRA respectively 

1-n = 4e^-2r?H + 2(l + c)e| + ^(3-5c)eH??H-^(3-c)C| + --- ; (6.1) 

= 6e^ - 2r]y - ^(44 - I8c)el - (4c - U)ey7]y - ^4 - ^(13 - 3c)e^ 
3 3d 

+ •••; (6.2) 

riT = -2e„ - (3 + c)e| + (1 + c)e„riH + • • • ; (6.3) 

= -2ey - ^{8 + 6c)el + ^{1 + 3c)eyr]y + ■ ■ ■ , (6.4) 

where c = 4(ln2 + 7) with 7 being Euler's constant. Notice the factors in 1 — n change 
even at first-order, due to the different definitions of r] which have been used. Similarly, we 
reproduce the second-order result for the ratio R of tensor and scalar amplitudes |^ 

25 

R = —e„[l + 2c{eH-r]H) + ---] (6.5) 



■^By contrast, the 'intermediate solution' can also be employed as the slow-roll solution in the simple 
potential V ex. (with /3 and cj) both positive), where the attractor hypothesis can be applied, though the 
solution to which the expansion process tends would have to be found numerically again. 
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25 



1 + 2 ( c - ^ ) (2e^ - r?v) + • • • 



(6.6) 



though it should be noted that this is not a direct observable [^] . Unlike the relatively sim- 
ple dynamics which we have emphasised in this paper, no way of extending these expressions 
analytically to arbitrary order is known. 

Final note: As we were completing this paper we received a preprint by Lidsey and Waga 
p2| which also discusses the slow-roll approximation, although with a considerably different 
emphasis. 
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Appendix 

We provide here a list of expressions, deemed too cumbersome and obtrusive to be imposed 
upon the main body of text, but which could prove useful in certain applications. The first 
four are exact extensions of Eq. ( p. 14 ) to higher-order parameters. 



riv = [3-e^]-^(3eH + 3r/H-r/|-£|), (A.l) 

£2 = [3-e^]-^{27eH7]H + 9f^-9eH7]l-l2r]HCl-3^l + 3vUl + VHCTl), (A.2) 
4 = [3 -e«]~^(81e^77| + 108eH£| + 274 -54eHr/|-72eHr?«£|- 5477^4 

- 27£| - 94 + 9e^r/| + Ue^vUl + 18r/^^£i + 27ry|cTi + 6r/^4 
-3vUh-Wh<^1-v14), (A.3) 

4 = [3- Eh]"' (810e«r/«e| + 4056^4 + 8l4-810e«r/|£|-405e« 7/^4 

- 270^14 - 216r,„4 " 27(1 + 270e^7?|e| + ISSe^vWH + 270r?«e|a| 
+ 1627714 + 27v„Cl - 30e„vUl " 15e«??|4 " ^^vUWh " 48r?|r| 

- 97?|Ci + Wt^UWh + + vIC'h) ■ (A.4) 

It is also possible to express any parameter we choose as a first-order differential relation in 
terms of lower-order parameters. This was done in Eqs. (|2T^ ) and (|2T^ ) for rj; we do this 
here for £ and a 



= 2ey7jy-J'^^r,'^, (A.6) 
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' m 



PI 



TT 



' m 



PI 



IT 



(A.7) 
(A.8) 



These compressions allow us to express any PSR parameter of order n as a first-order 
differential relation involving HSR parameters of order not exceeding n, as was done in 



Eq.( 2.15 ) for rjY- We present the case for ^v-, although such a result may be derived for any 
of the higher-order parameters, 



£2 



+ {vh - 3) 



PI 



IT 



{AS 



Eqs. ( A.lD "( [Alg| ) are all exact. We now give some approximate formulae, inverting some of 
the above relations to yield expressions for HSR parameters in terms of PSR parameters. 
If necessary, these can be fitted to Fade or Canterbury approximants, using the methods 



outlined in Section |5^. We have already stated the result for (Eq. (|5.3|) ); here we give 
the higher-order parameters. 



Vv 



2 



46 2 



17 



1 



13 



itr - ^ey7]y + 3e^ - 20et + 26e^r?v - 7ev< - y e^^v 



al - 3ey7]l + ISe^r?^ 



154 



(A.IO) 

(A.ll) 
(A.12) 



Note that these inversions are only valid when the attractor condition, Eq. (2.8), holds. The 



second-order truncation of Eq. ( A.IO ) is compatible with the result presented in |19|. 
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Figure Captions 



Figure 1 

A comparison of different analytic approximations with the exact, numericahy generated, 
solution H[(j)) for a potential V{4>) oc 0^, near the end of inflation. The normalisation of H 
is arbitrary. The slow-roll approximation and its first and second-order corrected versions 
are shown, together with the [1/1] rational approximant introduced in Subsection 5^, which 
is also a second-order correction. The rational approximant performs the best, as is more 
clearly seen in Figure 2. 



Figure 2 

The same comparison as Figure 1, but this time showing the exact enifp) corresponding to 
each of the solutions. Recall that the end of inflation is at e// = 1. The pathological behaviour 
of the first-order corrected solution, for which inflation never ends in this potential, is 
clear; all other solutions have a satisfactory end to inflation, with the rational approximant 
providing the best overall approximation to the exact solution. 
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